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Experiments on one- dimensional small capacitance Josephson Junction ar- 
rays are described. The arrays have a junction capacitance that is much 
larger than the stray capacitance of the electrodes, which we argue is im- 
portant for observation of Coulomb blockade. The Josephson energy can 
be tuned in situ and an evolution from Josephson-like to Coulomb blockade 
behavior is observed. This evolution can be described as a superconduct- 
ing to insulating, quantum phase transition. In the Coulomb blockade state, 
hysteretic current-voltage characteristics are described by a dynamic model 
which is dual to the resistively shunted junction model of classical Josephson 
junctions. 

PACS numbers: 73.40Gk, 74.50-hr 



1. INTRODUCTION 

The Coulomb blockade of Cooper pair tunneling (CBCPT) is a remark- 
able phenomena, where a normally superconducting tunnel junction becomes 
insulating, behaving as a capacitor with a critical voltage for current flow. 
The CBCPT can be observed when the Josephson coupling energy, Ej is the 
same order of magnitude as the charging energy, Ec = e^/2C, where C is 
the junction capacitance. The temperature must be low, ksT < Ec, Ej, and 
most importantly, the junction must experience a high impedance electro- 
magnetic environment. The electrodynamic environment controls the quan- 
tum fluctuations of the Josephson phase, and thereby determines whether 
superconducting (Josephson-like) or insulating (Coulomb blockade) behavior 
is observed. 

When the electrodynamic environment is described by the impedance 
Ze{u!), quantum fluctuations of the phase become large if Re[Ze{ui)] ^ Rq, 
where Rq = h/Ae^ = Q.ASkVt, is the quantum resistance for Cooper pairs. 
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Large quantum fluctuations of the phase mean that the number difference of 
Cooper pairs across the junction becomes a sharply defined quantum vari- 
able, as ANA(f) =1/2 for a coherent state. By realizing this high impedance 
environment with Josephson tunneling, we realize a unique state of charged 
matter, where many Bosons (Cooper pairs) are condensed into a state with 
well specified number. 

For a single Josephson junction in an arbitrary, linear electrodynamic 
environment, the effect of quantum fluctuations on Josephson tunneling is 
well described theoretically for the case Ej <ti Ec 0- However, when a single 
Josephson junction is placed in the environment of several other Josephson 
junctions, the theory is much less clear. Theoretical work on Josephson junc- 
tions arrays has primarily concentrated on two-dimensional (2D) arrays. A 
phase diagram is often calculated which maps out insulating and supercon- 
ducting regions, depending on the various parameters of the junctions, Ej, 
Eq, the normal state tunneling resistance R^, or the quasiparticle tunneling 
resistance, Rqp (for a review, see i). The theory can be related to experimen- 
tal work on granular thin films of superconducting metals (for a review see 
i) and 2p Josephson Junction Arrays fabricated with electron-beam lithog- 
raphy uB. Finite size effects, and the interplay between the island stray 
capacitance and the junction capacitance are generally not considered in the 
theoretical treatment of Josephson junction arrays. However, experiments 
on 2D arrays show a connection between the axray size and the character 
of the superconductor-insulator (S-I) transition Q. The Delft group has also 
studied long, two-dimensional arrays, which approach a ID parallel array 13, 
as well as short ID series arrays 

In this work we will examine experimental results on one-dimensional 
(ID) arrays of small capacitance Josephson junctions. We will describe the 
evolution from superconducting to insulating behavior in the arrays, which 
qualitatively is very similar to that observed in granular thin films and 2D 
Josephson junction arrays. We argue that the ratio of the junction capac- 
itance to the island stray capacitance is an important parameter for the 
coupling to a dissipative environment, and therefore is crucial for observa- 
tion of the S-I transition. 



1.1. One Dimensional Arrays 

Theoretical work on ID arrays of small capacitance Josephson junctions 
has analvzed the charge dynamics of ID arrays in the Coulomb blockade 
state liil Other theoretical analysis has centered on the questions of 

determining a phase diagram, and for what values of the array parameters an 
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insulating phase (Coulom 
effect) can be expected ^ 



blockade) or Superconducting phase (Josephson 
|i5l|i|,0 Theoretical analysis of small capacitance 
ID arrays usually begins with simplification of the charging energy. In many 
cases one takes only "on site" charging, where only the stray capacitance of 
each electrode to ground, Cq is considered. A better approximation to the 
experiments is to reduce the capacitance matrix to a bi-diagonal form, where 
Cq and the junction capacitance, C are taken into account. In experiments, 
typically C > Cq. In our experiments, we strive for C Cq by making 
the junction capacitance fairly large ( 1 to 3 fF) and the spacing between 
junctions very small (0.2/im) which leads to a smaller Cq. 

When C > Co the electrostatic screening length extends over several 
junctions. The potential due to one excess charge decays exponentially with 
the characteristic length A = ^C/Cq. This potential profile is known as a 
charge soliton II2 because the profile moves with the tunneling charge. Re- 
ducing the full capacitance to only C and Cq is valid only when a ground 
plane is located closer to the array than the junction spacing. One sample 
(B6-22 series described below) had a ground plane located 1.5//m from the 
array, and other samples had no ground plane. Thus the bi-diagonal capac- 
itance matrix is only approximate. In the absence of a ground plane, the 
potential distribution will decay less steeply than exponential, and go as 1/r 
at large distance. The charge soliton actually spreads out, and a slightly 
more effective polarization of the array occurs over the same characteristic 
distance, v'C/Cq 

Not only the electrostatics,.^ but also the electrodynamics of the array 
is effected by a large ^C/Cq ^3. For large ^C/Cq the junctions become 
strongly coupled, and the dynamics is dominated by Josephson plasmons. 
To demonstrate this point we consider an infinite ID array of Josephson 
junctions. For small phase gradient, when the current is less than the critical 
current, / <^ Ej, the Josephson phase-current relation can be linearized, 
and each junction is an effective inductance, Lj = ^q/IttIc, where <I>o = 
h/2e is the flux quantum. Replacing each Josephson junction by a parallel 
combination of Lj and C, we may calculate the impedance of the ID network 
shown in fig |^. The results of such a calculation are shown in fig. |^ for an 
infinite array, and a finite array with = 50 junctions, for typical values of 
Ic, C and Co in our experiments. 

For a finite array, we see a set of discrete resonances which are the 
Josephson plasmon modes. The width of these resonances is determined 
by the terminating impedance, Zq. For the infinite array, these resonances 
form a dense set of modes, where the impedance is much larger than Zq. At 
frequencies much less than the Josephson plasma frequency, ujp = \jSEjEc 
and for C 3> Co the real impedance can be written as. 
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Fig. 1. a) A schematic of the small capacitance Josephson junction array, 
b) An equivalent lumped element network when the phase difference across 
each junction is small. 




'^a{i^) = \j— = RQ\j—^\j— > iO<ujp = 

To observe the CBCPT, it is necessary that Ec ^ Ej and both Ec, Ej ^ 
ksT. If Ec 3> Ej the maximum Cooper pair current ("Zener current" ) is 
suppressed exponentially and becomes immeasurable, and if Ec <C Ej the 
threshold voltage is exponentially suppressed and becomes immeasurable 0. 
Thus, Za S> Rq requires that C ^ Cq. This impedance does not come from 
electromagnetic fields, as in a usual transmission line, but rather due to long 
wave-length Josephson plasmons involving the phases of many junctions. 

The large ratio y/C/Co ~ 10 in our arrays not only results in a large 
impedance, and therefore a Coulomb blockade, but also influences the effect 
of random offset charges. Simulations show that the static background po- 
tential due to random offset charge is smoothed by the averaging effects of 
the large ratio C/Cq @. One would expect that pinning of the charge in some 
local energy minimum in the array could be practically eliminated. However, 
simulations show that incoherent charge transport is substantially effected 
by static random background charge, even in 

the limit C> Co iliffl. To 

our knowledge, the effect of random offset charges on coherent Cooper pair 
tunneling has not been investigated theoretically. 
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Fig. 2. The impedance of the network shown in fig. |^ calculated in the limit 
C > Co for the parameters of sample B6-21 (see table I): C = 3.5x10"^^ F, 
Lj = 1.1x10-9, Co = 1.6xl0-^^ and ^Lq/Cq = 
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Fig. 3. Scanning Electron Microscope images of the arrays of SQUIDs show- 
ing how the array were terminated for the 2 point (upper) and 4 point (lower) 
measurement scheme. 
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2. EXPERIMENTAL RESULTS 

The arrays are actually serially connected SQUIDs (Superconducting 
Quantum Interference Devices). The SQUID geometry allows one to exter- 
nally impose a phase shift between neighboring superconducting electrodes 
by application of a magnetic flux. This external phase shift can be consid- 
ered as effectively tuning the Josephson coupling energy with an external 
magnetic field, Ej{B). Figure ^ shows scanning electron micrographs of 
two different arrays. The electrodes were Al, deposited on an oxidized Si 
substrate, with an AI2O3 tunnel barrier. At the right hand side of fig.|3| we 
see how the arrays were terminated. Two different arrangements were used. 
Figure ^a shows a 2-point termination, where a wide strip connects to the 
first SQUID. This 2-point arrangement had the disadvantage that expelled 
magnetic fiux from the wide strip focused more flux into the edge loops. Con- 
sequently the two edge junctions had a smaller period of modulation with 
magnetic field. Figure ^ shows the 4-point termination, where no edge ef- 
fect with the magnetic field was observed. For the 4-point termination, the 
junctions in the leads are not SQUIDs, and thus do not have a tunable Ej. 

Table |l] shows the main parameters of the arrays studied thus far. Sam- 
ples with the same chip code (B6-##) were made simultaneously on the 
same chip. The samples are arranged with decreasing Ejq/Ec- The note 
indicates whether the array was terminated with a 2 point or 4 point lead 
configuration. The value of -Ejo = {^q/2t[)Ic = Rq/S.q/2R]\j was determined 
from the measured normal state resistance Rn and the superconducting en- 
ergy gap, Aq. The charging energy Ec = e^/2C was calculated from the 
observed area of the junction and the specific capacitance Cg = 45fF//im^. 
Depending on the array parameters, different behavior of the arrays could 
be observed. We divide these into three different regions discussed below. 



2.1. Superconducting State 

For the samples in table || with a smaller ratio Lj/Cq, the current- 
voltage (I-V) characteristic exhibited Josephson-like behavior at zero mag- 
netic field. Figure ^ shows the I-V curve of sample B6-21a. Here we see a 
nearly vertical "super current" branch of the I-V curve near zero voltage, 
and a successive set of vertical quasi-particle branches separated in voltage 
by 2Ao/e. The dc load in the bias circuit was such that we could not trace 
out the negative differential resistance between these vertical quasi-particle 
branches. At higher voltages the successive quasi-particle branches merge 
into one flat branch with a current independent of voltage. At even higher 
voltages, the I-V curve eventually switches to the normal tunneling branch. 
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Table 1. Main parameters of studies arrays. The notes have the following 
meaning: 2 pt. or 4 pt. describes the array termination. Sl-trans means that 

superconductor-insulator transition could be observed. Short arrays showed 
no Coulomb blockade of Cooper pair tunneling (CBCPT) and longer arrays 
with smaller Ejq/Eq showed only CBCPT (no superconducting state) 



Sample 


N 


Rn 


Ejo 


Ec 


Ejo/ Ec 


VLj/Co 


Note 




code 


Jun. 


(kO) 




ifieY) 






(kO) 






B6-21a 


255 


1.1 


586 


23 


25 




8.5 


2 pt., 


Sl-trans 


B6-21b 


127 


1.1 


586 


23 


25 




8.5 


2 pt., 


Sl-trans 


B6-21C 


63 


1.1 


586 


23 


25 




8.5 


2 pt.. 


Sl-trans 


B6-21d 


31 


1.1 


586 


23 


25 




8.5 


2 pt.. 


no CBCPT 


B6-21e 


15 


1.1 


586 


23 


25 




8.5 


2 pt.. 


no CBCPT 


B6-21f 


7 


1.1 


586 


23 


25 




8.5 


2 pt.. 


no CBCPT 


B6-41a 


255 


1.6 


406 


29 


14 


10.2 


4 pt.. 


Sl-trans 


B6-32a 


255 


2.7 


227 


29 


7.8 


13.3 


4 pt. 


Sl-trans atRg 


B6-22a 


255 


4.6 


142 


23 


6.1 




17.4 


2 pt.. 


Sl-trans 


B6-22b 


127 


4.6 


142 


23 


6.1 




17.4 


2 pt.. 


Sl-trans 


B6-22C 


63 


4.6 


142 


23 


6.1 




17.4 


2 pt.. 


Sl-trans 


B6-41b 


255 


4.3 


149 


50 


3.0 


16.8 


4 pt.. 


only CBCPT 


B6-32b 


255 


8.5 


78 


52 


1.5 


24 


4 pt.. 


only CBCPT 


K9-32 


255 


10.4 


62 


59 


1.1 


26 


4 pt., 


only CBCPT 


K9-43a 


65 


15.4 


42 


59 


0.71 




32 


4 pt.. 


only CBCPT 


K9-43b 


255 


15.7 


41 


59 


0.70 




32 


4 pt.. 


only CBCPT 
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Fig. 4. The current voltage characteristic for sample B6-21a measured at 
T = 50mK and B = 0. 



The I-V curve shown in fig. ^ is not that expected for arrays of classical, 
non-interacting Josephson junctions. As the bias is increased in such arrays 
we would expect an increasing series of switching currents, starting from the 
lowest critical current in the array. For these small capacitance, strongly 
coupled Josephson junctions, we rather observe a decrease of the switching 
current as we move out in voltage, as well as an increase of the retrapping 
current. We have no quantitative explanation for this observed behavior. We 
speculate that it has to do with electromagnetic interaction of the junctions, 
which are strongly coupled due to the very small Cq. 

As the magnetic field is increased, the critical current of all junctions 
is decreased. The slope of the "critical current" branch also increases, and 
as Ej is further suppressed, a distinct threshold voltage emerges as the 
Coulomb blockade of Cooper pair tunneling (CBCPT) becomes observable 
(see fig. 1) 



2.2. Superconductor - Insulator Transition 

Figure ^ shows three I-V curves at different magnetic fields. The nearly 
vertical line (I-V at B = 57G) has a finite slope corresponding to about 
300 kJ7. The critical current and the hysteretic switching between the "su- 
percurrent branch" and the quasiparticle branch is off scale for this curve. 
As the magnetic field is increased the I-V curve develops a very high resis- 
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Fig. 5. The I-V characteristics at T = 50mK for three different magnetic 
fields as the array B6-22a goes through the Superconductor-Insulator tran- 
sition. The curve labeled 57G is actually on the insulating side of the tran- 
sition, as Rq increases when T — > at this magnetic field. 



tance state for voltages below a threshold voltage, which is the CBCPT. The 
curve labeled B = 70G in fig ^ shows a hysteretic I-V curve in this Coulomb 
blockade state, meaning that "back bending" or region of negative differen- 
tial resistance exists at low current. We will discuss this hysteresis in the 
next section. The intermediate curve {B = 66G)in fig. |5| shows remnants of 
both types of hysteresis. 

Figure |5| shows that the crossover to the CBCPT is characterized by 
a set of very nonlinear I-V curves. We would like to know the linear re- 
sponse of the array as in undergoes this transition from Josephson-like to 
Coulomb blockade behavior. To this end we have devised a measurement 
of the zero-bias resistance, i?o = limv dVrms/dlrms under the condition 
that the power dissipation by the array is constant. Rq is measured by 
phase sensitive detection of a small signal excitation with two lock-in am- 
plifiers. The excitation is adjusted so that the product of the two signals 

{dI)rn,s{dV)rms = lO'^^ Watt. 

Rq could be measured versus temperature and magnetic field. Figure ^ 
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displays the results of such measurements. At low magnetic fields (bottom 
curves of fig. Rq decreases as the array is cooled below the transition 
temperature of the Al (Tc = 1.2K) and eventually becomes flat, showing 
a temperature independent resistive state as T ^ 0. The finite resistance 
corresponds to the slope of the " critical current" branch of the I-V curve (see 
curve B = 57 in fig. |5|). Experiments indicate c3 that in longer arrays the 
resistance of this "flat tail" decreases, indicating the non-zero resistance in 
this "superconducting" state is due to the finite size of the array. When the 
magnetic field is increased, the resistance rises and the "flat tail" remains. 
Further increasing the magnetic field causes a transition to a different kind 
of behavior, where Rq increases as T — > 0, indicating the development of a 
Coulomb blockade. 

We have previously shown IH3 how length scaling of Rq could be used 
to determine the critical point between superconducting and insulating be- 
havior. Length scaling analysis allowed us to determine a T = critical 
point where Rq is independent of length. The critical point separates those 
curves with a "fiat tail" in Ro{T) from those with increasing Rq as T ^ 0. 
In the experiment shown in fig. |6| we find that this transition occurs at 
Rq{T = 0) ~ h/Ae^ = 6.45kr2 (see the arrow in fig. P). The transition 
does not always occur at Rq , and the measurements displayed in fig. ^ were 
special in that the array was terminated with a 4 point configuration. The 4 
point termination should reduce any series measurement of dissipation aris- 
ing at the edge of the array. However, one junction is probably not enough 
to effectively isolate the array from the environment. Further experiments 
are needed to determine if this transition at Rq = Rq is universal. 



2.3. Insulating State 

When the CBCPT is well developed in the arrays, Rq is essentially in- 
finite, and we have a zero current state below the threshold voltage. Figure 
shows how the threshold voltage changes with magnetic field on the in- 
sulating side of the S-I transition. We typically find that Vt{B) exhibits a 
sharp cusp as Ej when ^ext — *■ ^of^ {B = 78 G for this sample). 
For samples B6-22 and B6-21 where we could measure length effects, we 
found that Vt was roughly proportional to length. The threshold voltage de- 
pends on the magnetic field in a periodic way, with period corresponding the 
fiux quantum in each loop. The magnitude of the current is also periodic 
in the magnetic field. These observations prove that the current and the 
threshold voltage are due to Cooper pair tunneling. Although the Cooper 
pair tunneling is a coherent process, we measure dissipation (finite current 
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Fig. 6. The S-I transition is seen by measurement of the zero bias resistance 
vs. temperature. The data are taken on sample B6-32a at several values of 
the magnetic field in the range to 64 G. 
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Fig. 7. The magnetic field dependence of the threshold voltage for sample 
B6-22. The inset shows the I-V curves. 



and non-zero voltage) because the pair tunneling interacts with dissipative 
degrees of freedom. 

As in the case of a single small capacitance Josephson junction, dissi- 
pation arises fromexcitation of the electromagnetic environment or quasi- 
particle tunneling 0. Quasi- particle tunneling becomes important when the 
voltage drop across any junction exceeds the gap voltage, 2A/e = 400fjM. 
It is important to note that although the voltage over the entire array may 
be greater than 2A/e, this voltage is distributed over many junctions. The 
charge soliton length for Cooper pair charge solitons is determined by an 
effective capacitance of the junction 0, As = Cg/ //Co = \/{2e/2TTVc)/Co. 

The critical voltage Vc{Ej / Ec) goes to zero exponentially for Ej ^> Ec 0. 
Hence, — > oo at the S-I transition and near this transition the voltage is 
dropped uniformly over the entire array. 

In the insulating state, a very interesting model for charge transport 
has been developed, which has a direct duality to flux transport in long 
Josephson junctions 0. A static version (/ = 0) of this model has been used 
to explain the magnetic field dependence of the threshold voltage 113. We 
have been able to extend the model to the time domain, calculating the dc 
I-V curve and fitting it to the data. This dynamic model is able to explain 
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the hysteresis observed in the dc IV curve (see fig. |8| and the curve labeled 
B = 70G in fig. 

In the insulating state, the Josephson junctions in the ID chain are de- 
scribed in terms of the dimensionless quasi-charge x- If ^ phenomenological 
resistance R and inductance L are introduced together with the capacitance 
Co, a sine-Gordon like model can be derived to describe the dynamics of 
x{x, Soliton solutions exist for this model, and the inductive term is 
what gives the charge soliton mass, or inertia. We may simplify the full 
sine- Gordon-like model by dropping the second spatial derivative in order 
to capture the basics of the dynamic behavior. The array is then viewed as 
one lumped element, which is valid if the array is shorter than the soliton 
length, as is the case near the S-I transition. Alternatively, the second time 
and spatial derivatives may combined by assuming a traveling wave solu- 
tion to the sine-Gordon-like model. Either approach results in the following 
equation c3: 

V + RIz . , , 

77 = PX + X+ saw (x) 

Vc 

where V is the applied voltage, and the dc current is given by the time 
average < x>. Here saw(x) is a 27r periodic function describing the charging 
and discharging of the junctions capacitance by single Cooper pair tunneling 
events. The dot refers to differentiation with respect to dimensionless time, 
t/ RC. For the lumped element case, the dimensionless parameter. 



R^ 



2e 



2-kVc 



determines the amount of hysteresis observed in the dc I-V curve. We have 
modeled the dissipation arising from Zener transitions in the Bloch energy 
bands of the single junction Coulomb blockade Bj. To make the model simple 
we assumed that the resistance is zero below the Zener threshold current, 
Iz, and characterized by a linear resistance R above Iz- 

The parameters Vc, R, and Iz are fixed for a particular magnetic field 
(Ej/Ec), and are determined from the observed I-V curve. The parameter 
/3 or L is then adjusted to fit the observed amount of hysteresis. Good fits 
to the data can be made, as seen in fig. ^. However, the inductance L 
required to explain the observed hysteresis is extremely large, ~ 10~^ H. 
It is certainly not due to electromagnetic inductance (~ 10~^^ H) and too 
small to be due to Josephson inductance {N^o/2itIc ~ 10~^ H). We do 
observe that the value of L obtained fromthese fits does diverge as /c — > 
as expected for the Josephson inductance c3. 
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Fig. 8. Hysteretic I-V curves in the insulating state are observed. The 

open circles are fits to the experimental data (dots) using the non-linear 
dynamic model for the quasicharge which is dual to the RS J model of typical 
Josephson junctions. 
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Although the origin of this large inductance is not understood, we em- 
phasize that an inductive term is absolutely necessary if we are to explain the 
hysteresis with any dynamic model. We further stress that this hysteresis 
is clearly associated with the presence of Josephson Coupling in the arrays. 
Such hysteresis has not been observed in any normal tunnel junction arrays 
to the best of our knowledge. In our arrays, the normal state behavior shows 
no clear zero current state. The normal sate is achieved at large magnetic 
fields {B ~ lOOOG) where the superconducting energy gap is suppressed. 
Only a weak remnant of the Coulomb blockade exists in the normal state 
because the normal tunnel junction resistance is the order of the quantum 
resistance. 



3. CONCLUSION 

We have given an overview of the measured transport characteristics 
of ID small capacitance Josephson junction arrays. The interplay between 
Josephson and Charging energies is clearly seen as a transition between 
Josephson-like to Coulomb blockade I-V curves. Coulomb blockade is ob- 
served even when Ej > Ec, depending on the length of the array. We 
point to the importance of the stray capacitance Cq and the long wave- 
length Josephson plasmons in determining the coupling to dissipation. A 
quantum phase transition can be observed by measuring the temperature 
dependence of the zero bias resistance as the effective Josephson coupling 
is modulated in the arrays. In the insulating state, the hysteretic current 
voltage characteristics can be explained by a dynamic model which is dual 
to the usual RSJ model for Josephson junctions. 
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